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We study the dynamics of open quantum many-body systems driven across a critical point by 
quenching an Hamiltonian parameter at a certain velocity. General scaling laws are derived for the 
density of excitations and energy produced during the quench as a function of quench velocity and 
bath temperature. The scaling laws and their regimes of validity are verified for the XY spin chain 
locally coupled to bosonic baths. A detailed derivation and analysis of the kinetic equation of the 
problem is presented. 



I. INTRODUCTION 

A series of beautiful experiments on the dynamics 
of cold atomic gasesi*^!^ spurred renewed interest in 
the study of non-equilibrium quantum many-body sys- 
tems. On the theoretical side these experiments triggered 
an intense investigation mostly devoted to the simplest 
paradigm of nonequilibrium quantum dynamics: the con- 
trolled variation in time of one of the system parameters 
(quantum quenches). In the case of sudden quenches, 
where the driving parameter is changed on a time scale 
much shorter than typical time scales of the system, a 
number of important issues have been addressed. We 
mention, for example, the study of the signatures of 
universality in the quench dynamics of quantum criti- 
cal systems^, the presence of thermalization in integrable 
vs. nonintegrable systems^, as well as the description of 
generic nonequilibrium quenches using thermodynamic 
variables — and their statistics^. 

In this paper we will focus on the opposite case in 
which the control parameter is varied adiabatically, a 
case which becomes particularly interesting when a crit- 
ical point is crossed during the adiabatic evolution. Be- 
cause of the vanishing of the energy gap at criticality, 
the system is unable to follow adiabatically the driving 
remaining in its equilibrium/ground state when passing 
through the quantum critical point the system will not 



be able no matter how slow is the quench. The study of 
these deviation from the adiabatic dynamics is a prob- 
lem which is very important in a number of different 
branches of physics ranging from the defect formation 
in the early universe^i^ to adiabatic quantum computa- 
tion — or quantum annealingiiii^. Depending on the 
context, the loss of adiabaticity has been characterized 
by the excess energy at the end of the quench, by the 
density of defects (if the final state was a fully ordered 
system), or by the fidelity of the time evolved state with 
the ground state of the Hamiltonian at the end of the 
quench. 

The scaling of the density of excitations gener- 
ated during the dynamics as a function of the ve- 
locity of the quench was first predicted in Ref. [Tsl . 
[14] for a quantum critical system. The mecha- 
nism behind the generation of excitations/defects is 
similar to so-called Kibble-Zurek (KZ) mechanism^ 
first proposed for classical phase transitions. Fol- 
lowing these initial works a number of specific mod- 
els were scrutinizedi^^^^'ia^'21^^^^^^^^, 
thereby confirming the general picture. 

All the works mentioned previously assumed unitary 
Hamiltonian dynamics. We know, however, that under- 
standing the effect of the external environment on the 
adiabatic dynamics is of paramount importance for sev- 
eral reasons. In the case of adiabatic quantum computa- 
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tion, decohcrence is a fundamental limiting factor to the 
ability of implementing quantum algorithms. Further- 
more, an experimental verification of the KZ scaling in 
a quantum phase transition can only occur through the 
detection of this effect at low temperatures, i.e. when 
the quantum critical system is in contact with a thermal 
bath. Despite its importance the adiabatic dynamics of 
open critical systems is a much less studied problem. The 
effect of classical and quantum noise acting uniformly 
on a quantum Ising chain was considered in Reflsol and 
ReflSlI respectively. Numerical simulations for a model of 
local noise on a disordered Ising model were performed 
in Refill. Moreover the effect a static spin bath locally 
coupled to an ordered Ising model is studied in ReflSSl. 
In a recent Letter—, we have addressed the universality 
of the production of defects in the adiabatic dynamics in 
the presence of an environment by generalizing the scal- 
ing theory to open critical system and by formulating 
a quantum kinetic equation approach for the adiabatic 
dynamics across the quantum critical region. We found 
that, at weak coupling and for not too slow quenches the 
density of excitations is universal also in the presence of 
an external bath. In this paper we extend the results 
presented in Ref. |33 and provide a detailed derivation of 
the kinetic equations and of the scaling approach. 

The paper is organized as follows. We first derive qual- 
itatively the scaling laws obeyed by both the density of 
defects and of energy generated in a quench for a generic 
open quantum critical system (Sec|TT|. We then address 
a specific one-dimensional model possessing a quantum 
critical point: the XY spin chain in transverse magnetic 
field. To model a thermal reservoir we couple the system 
to a set of bosonic degrees of freedom, as in the spin- 
boson model. Baths are chosen with power-law spectral 
density and are locally coupled to strings of neighboring 
spins. The model, a generalization of the one studied 
in Refls^. is discussed in Scc lIIII For this model, we 
derive a kinetic equation within the Keldysh technique 
(Sec. IIVI and Appendix A1-A2) which allows us to com- 
pute the density of defects. In Sec. |V] (and Appendix B) 
we discuss the spectrum of relaxation times needed for 
a comparison with the scaling approach. The density of 
defects and of excitation generated in a quench, and a 
comparison with the scaling laws is presented in Sec. IVII 
Finally in Sec lVIII wc summarize our conclusions. 



II. SCALING ANALYSIS 

In this section wc discuss the scaling laws obeyed by 
the density of excitations^ and by the energy density fol- 
lowing a linear quench of a control parameter h from an 
initial value hi to a final one hf, through a second-order 
quantum critical point at he- The system, during the 
whole dynamics, is kept in contact with a bath at tem- 
perature T. In the h — T plane the adiabatic quench is 
described by the horizontal line shown in Fig[T] For adi- 
abatic quenches occurring at zero temperature, the sys- 



tem stops following the external drive adiabatically and 
can be considered as frozen around the quantum criti- 
cal point. This happens roughly when the time it takes 
to reach and cross the quantum critical point becomes 
comparable to the internal time scale (the inverse gap 
A(/i) ~ \h — hc]"'^^). The determination of this crossover 
point is the fundamental ingredient which leads to the 
scaling in the case of unitary cvolutio n^^i^^ . In the case 
of a finite temperature quench there is a new important 
timescale which enters the problem, the time at which 
the system enters (and eventually leaves) the quantum 
critical region (see FiglT|) . Initially the system is in equi- 
librium with the bath at a low temperature (T ^ A(/ii)) 
and the behavior is, to a large extent, as in the zero- 
temperature case. In the quantum critical region^, char- 
acterized by the crossover temperature T ~ |/i — /id'^^, 
the gap is much smaller than the temperature itself. One 
can therefore expect that during the interval the system 
spends in the quantum critical region a number of exci- 
tations will be produced by the presence of the environ- 
ment. Interestingly also this contribution to the defect 
production obeys a scaling law^. 
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Figure 1: A sketch of the finite temperature crossover phase- 
diagram close to the quantum critical point. Crossover lines 
T ~ \h — hc\"^ separating the semiclassical regions from the 
quantum critical region are shown. The latter is traversed by 
the system during the quench in a time Iqc ■ 

We now proceed with the derivation of the scaling laws. 
In the rest of the paper we will consider the density £ and 
the energy density E of excitations, defined respectively: 



S = 



E 



(2^ 



EuVk 



(1) 
(2) 



where Vk is the population of the excitation with quan- 
tum number k, Ek the energy spectrum at hf, and d is 
the dimensionality of the quantum system. 

The first assumption we make consists in separating 
the density of excitations/energy at the end of the quench 
in the sum of two (coherent and incoherent) contribu- 
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tions: 

^ — ^coh ~1~ ^inc : (3) 
E ~ Ecoh + Einc ■ (4) 

In the previous equation. Ecoh i^coh) is the density of 
energy (excitations) of the system produced coherently 
in the absence of the bath; the incoherent contribution 
Einc {£inc) arises instead from the bath/system interac- 
tion. The separation of a coherent and an incoherent 
contribution, eqs. ([3]) and (jl]), requires weak couphng a 
between the system and the bath. Evidence for the vahd- 
ity of this assumption wiU be shown below (see Eqs.([TU|) 
and 111])). 

In the absence of an environment, the density of exci- 
tations was shown to obey the KZ scalin g^^i^^ 

£coh = £kz cc . (5) 

In order to obtain a similar relation for the energy density 
in Eq. (|30p additional information on Ek at hf is needed. 
Thus, the scaling of this quantity depends on the details 
of the system at the end of the quench. A simple scaling 
law can be obtained only in specific situations, e.g. for 
quenches halted at the critical point hf = he, where Ek oc 
fc^. By using techniques similar to those employed in 
Ref. [ij one obtains 

Ecoh oc . (6) 

Let us now derive a scaling law for the incoherent con- 
tributions Sine and Einc- To this end it is convenient 
to divide the quench in three steps (see Fig. [J): ini- 
tially the system is in the so-called low-temperature re- 
gion at T <C A. Here the relatively high energy gap 
suppresses thermal excitation and the system remains in 
the ground state. Close to the critical point the system 
passes through the quantum critical region: thermal exci- 
tations are unavoidably created because of the relatively 
high temperature T A. As we shall see below, the 
density of excitations generated in this region arc uni- 
versal functions on the velocity of the quench and on 
the temperature as long as only the low-energy details of 
the system spectrum matter. On the contrary, the bath- 
induced relaxation occurring once the system leaves the 
quantum critical region, entering the other semiclassical 
region (T <C A), depends on the details of the energy 
spectrum; hence the relaxation towards an asymptotic 
thermal state at temperature T is not expected to be 
universal if the final hf is far off the critical point he- In 
our analysis below we will neglect the effects of this non- 
universal relaxation. We are therefore assuming that the 
time elapsed between the moment when the critical re- 
gion is left and when the quench is stopped (and the mea- 
surement of the density of excitations/energy is made) is 
short as compared to the typical relaxation times in the 
semiclassical region. In Section IVII wc will further com- 
ment on this point for the specific case of the quantum 
XY chain, showing that the scenario just depicted holds 



for a wide range oi hf and v. The dynamics of the prob- 
ability to excite the model k Vk can be described, inside 
the quantum critical region, in terms of a phenomenolog- 
ical rate equation: 

j^n = -^[Vk-Vl'^{hc)] , (7) 

where V^{hc) is the critical thermal equilibrium distri- 
bution and is the relaxation rate, oc oT^ . As 
shown in Section Ivl and appendix C, 9 can be related to 
characteristics of the bath and to the critical indices of 
the phase transition (see Eq. (|56|) '). From the relation 
T ~ A ~ \h — hc\^^ we deduce the time spent inside the 
quantum critical region is 

A direct integration of Eq. ([7]) gives for the thermal ex- 
citation created in the quantum critical region Vk ^ 
(1 — e*Q^/'^)7^^'' (he). Finally, integrating the latter over 
all k-modes we get: 

S.nc OC (l-e*"^/^) jdEE^'^-^Vl^ihc) , (8) 

where we used the scaling of the excitation energy E oc 
. For the density of energy, a similar relation holds in 
the case of quenches halted at the critical point 

E,nc OC (l - jdEE^'^Vfihc) , (9) 

where tqc/'i is due to the fact that in this case only half 
of the quantum critical region is crossed. Finally, since 
the thermal distribution is a function of E /T , a simple 
change of variable gives the required results 

£,nc OC av-^T<'+^ , (10) 

E,nc OC av-^T'>+^^ , (11) 

that are vahd in the limit T'^/"^ < vt. Eqs. HO]) and 
PT|) together with Eqs. (O and © give, through the 
assumptions Eqs. ([3]) and (j4|), the general scaling-law for 
the quench dynamics of open systems. The different scal- 
ing of the two contributions with respect to the velocity 
V implies that for slow quenches v < Vcross the incoher- 
ent mechanism of excitation dominates over the coherent 
one, and viceversa for v > Vcross- The crossover velocity 
can be deduced by equating £^nc — £coh and Emc — Ecoh 
yielding: 

v^cross OC T^^+Z^+lni)i^+^) , (12) 

v^^^,, cx a -(^^ t(^+-("-"wi)(i+77i) . (13) 

III. QUANTUM XY MODEL WITH THERMAL 
RESERVOIR 

The scaling laws derived above will be tested against a 
specific model: an XY-chain coupled to a set of bosonic 
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baths. The Hamiltonian of the XY chain is defined as 

N / 

i ^ ^ 

(14) 

Here N is the number of sites (cr^' ^' ^ are Pauh ma- 
trices). Each spin is coupled to its neighbors d by 
anisotropic Ising-like interaction and subject to a trans- 
verse magnetic field h (the coupfings are expressed in 
terms of the exchange energy). In the thermodynamic 
limit ^ oo, a quantum phase transition at ft-c = 1 
separates a paramagnetic phase for h > 1 from a ferro- 
magnetic phase (/i < 1) where the Z2 symmetry is spon- 
taneously broken and a magnetic order along x appears, 

The spin Hamiltonian Eq.i \1A^ can be diagonalized 
by using the Jordan- Wigner transformationfiZ to map 
the spins into spinless fermions Cj and thus obtain in 
momentum space (after a projection in a definite parity 
subspace) 

fc>0 

Tik ~ — (coa k + h) Tz + J sink Ty , (15) 

where 'ij, = ( c[. C-k) ai'e Nambu spinors and f are 
Pauli matrices in Nambu space. Finally, a Bogoli- 
ubov rotation diagonalizes the Hamiltonian: Hs = 

Efe>o Afc(^fc'?fe - V-kvlk)^ where 

Afc = \/(cos k + h)^ + (7 sin fc)2 , (16) 

is the quasi-particle dispersion. At h = he the spectrum 
becomes gapless with a linear dispersion relation cx 
TT — fc; accordingly, the critical indexes of the model are 
V = z = 1. 

The spins are also locally coupled to a set of N/l 
bosonic baths 

N/l-l n-i \ 
j=Q \r=0 / 

where Xj = J^p ^isi^pj + ^-/3j) and bj^ j (6^j) are the 
creation (annihilation) operators of the j-th bosonic bath. 
As a result of the coupling in Eq. P?|) each baths corre- 
lated in a string of / adjacent spins. The model presented 
above generalizes the one considered in Ref. [13, where the 
spins were individually coupled {I = 1) to Ohmic baths 
(. = 1). 

The total Hamiltonian reads: 

H = Hs + H,nt + Hb (18) 

where Hb = p^/ib^pj^iij- The spectral density of the 
baths J(uj) — X]/3 A^(5(a; — LOp) is 

J(w) = 2auj'e-'^''^^e{uj) (19) 




Figure 2: A cartoon of the spins-baths coupling (|17|) for the 
I — 3, where each bath is coupled to three spins. 



where a is the system/bath coupling, ujc is a high energy 
cutoff and 9{s) is the step function^. 

In the situation we are interested, the system is initial- 
ized in its ground state at large h. The coupling of the 
spin to the bath through preserves parity symmetry 
(see Eq. (|20|) ). Therefore, once the system has initially a 
specific parity, it will remain in the corresponding sector 
for the entire evolution. Throughout this paper we con- 
sider N even: in this case the ground state has always 
an even number of fermions Ck and we are thus allowed 
to select the even parity sector and neglect the odd one. 

In momentum space bp^q = 

ES"^ exp(-ig.7) with g = ^ and af- 
ter the Jordan- Wigner transformation we get 

^™* = -;4EE^(9)^If^*fc+f^, (20) 

k q 

where F{q) = l/\/I^^~pCxp(— irg//). For I ~ 1 each 
spin interact with a different bath and according to 
Eq. (PO]) all k modes are coupled (i.e., transitions k ^ k', 
yk,k' are induced). In the opposite case of I = N (just 
one bath for the whole system) no transition between dif- 
ferent k is allowed. In the intermediate case each mode 
interacts with the other modes in an interval of width 
2Trl/N. 

It is important to notice that correlations between 
baths over a finite distance would not change qualita- 
tively our picture as long as we focus on the critical 
properties of the model. Indeed, near criticality the di- 
vergence of the correlation length makes the details of 
the bath correlations over microscopic distances unim- 
portant. For the same reason, as long as one is interested 
in the low-T properties of the bath, the specific value of I 
is not relevant provided l/N — > in the thermodynamic 
limit. Specifically, for all values of / such that T ^ ( j) 
the same dissipative dynamics is obtained, since transi- 
tions with large Ak (and hence large energy) are ther- 
mally suppressed (we used cx A:^ at a fixed T). In 
this regime, therefore, the system cannot resolve the mi- 
croscopic details of different system-bath couplings (i.e., 
whether / = 1,2,3,...). In the following, we thus focus 
on the case Z = 1: specific high-temperature and non- 
critical behaviors for different I could be easily investi- 
gated within the same scheme considered below. Only if 
I = N the dynamics of the system changes qualitatively. 
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IV. KINETIC EQUATION 

111 this section we derive a kinetic equation for the 
Green's function of the Jordan- Wigner fermions within 
the Keldysh formaUsm. In terms of this Green's func- 
tion we will then calculate both the excitation and the 
energy densities, Eq. ^ and Our analysis in terms 
of a kinetic equation will provide support for the scaling 
laws obtained above, Eqs. P^ - PT|) . while allowing us 
to study the non-universal dynamics beyond the limit of 
applicability of the scaling approach. 

The fermionic Keldysh Green's function is a matrix in 
Nambu space defined by 



spectrum of the systero^^. As previously mentioned, we 
will not consider this case here. Therefore, we assume 
that the bosons have a non-zero inverse lifetime F ^ T 
which provides a natural cutoff-time for the bath correla- 
tion functions. Within this assumption it is now possible 
to describe the kinetics of the system using a Markov ap- 
proximation together with a self-consistent Born approx- 
imation. The latter is justified for weak system/bath 
coupling {a <^ 1) and is represented diagrammatically 
in Fig. [3]- (a). We will neglect the tadpole diagram (b), 



(21) 



(see appendix Al for notations) where 7 is the Keldysh 
contour. In the following we will neglect the initial corre- 
lations between system and bath^. Hence 7 consists of 
just a forward and backward branch on the real time axis. 
Below we will sketch of the main steps of the derivation: 
the remaining details can be found in Appendix A1-A2. 

The starting point of our derivation is the Dyson's 
equation in its integro-differential form 



idt.-Hkiti) Guituh) = Sih-h) (22) 

dfSfe(ti,i)Gfc(f,t2) 



and an analogous one obtained by differentiation with 
respect t2- Here Sfc is the self-energy associated to the 
interaction of the system with the bath. In order to com- 
pute the energy and excitations density we need to find 
the equal-time statistical Green's functions. The latter 
are defined as 

[GtihMh,, = *(*l,(i2)*,.,(ti)) (23) 



[Gl{t,M)],a = -*(*fe.(ii)*L(i2)) 



(24) 



An equation for these correlators can be obtained from 
Eq. (|22|) by using standard techniques^ (see Appendix 
Al for details). For the equal-time Green's function 
G^{t, t) we obtain: 



idtG< = 



Hk^Gl: 



I]> • G< - S< • G> + G< • S> 



(25) 



Gk '^k ■• 



where the dots indicate the convolution: 



dtj:>{t,t)G<{t,t). 



In order to proceed with the solution of Eq. (^5)1 it is 
now important to discuss the approximations we make for 
the self-energy. Let us first notice that long-time corre- 
lations induced by the bath may change the universality 
class of the transition, by renormalizing the low energy 



a) 



b) 



Figure 3: Lowest-order diagrams contributing to the self- 
consistent Born approximation: dashed lines correspond to 
the non-interacting bath Green's function g, while solid lines 
to the interacting system Green function G. a) corresponds 
to Eq. (gSl), while b) to Eq. (|47)l . 

which represents just a small shift of the energy levels. 
By going to the interaction picture 



Gfc(ti, <2 



ul(h)Gkih,t2)Uk{t2) , 



it is now evident that, within our assumptions, the evolu- 
tion of Gfc can be considered slow as compared to that of 
the bath correlators appearing in the self-energies. Us- 
ing this separation of time scales it is possible to im- 
plement the Markov approximation and transform the 
general integro-differential kinetic equation into a simple 
differential equation (see Appendix A1-A2). We then ob- 
tain, in the case in which each spin is coupled to its own 
bath (/ = 1), the kinetic equation 



dtG< 



lJ2r^l + ^G<)D,kGl 



N 
where 



I [Hk. G<] = (26) 
T''G<bl^(l + iG<) + H.c. 



qk=i f dsg>{s)Ul{t,t- s)rUk{t,t~ s) , (27) 




g^{t) = —i{Xq{t)Xq{0)), and Uk{to,t) is the evolution 
operator satisfying idtUk — 'Hk{t)Uk- The left-hand-side 
of Eq. (p6|) represents the free evolution term, while the 
right-hand-side describes the scattering between the k 
modes mediated by the bath degrees of freedom. Notice 
that the number of equations scales linearly with the sys- 
tem size N, in contrast to conventional systems of master 
equations whose number scales exponentially with N as 
a result of the fact that the full density matrix (i.e. all 
TO— points Green functions) is considered. The fact that 
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we are considering only the two-point Green's function 
self consistently using the Born approximation is respon- 
sible for the non- linear nature of Eq. ([26]) , in contrast to 
the linearity of the master equation. 

In the eigenbasis of the Hamiltonian Ti^ the Green's 
function can be parameterized as 

where Vk — {Vk^k) is the population of the excited mode 
k and Ck = (rj-kVk) can be regarded as a "coherence" 
term^. In the static case, where the evolution operator 
isUk = exp{—i'Hkt), the stationary solution of the kinetic 
equation (|26p is correctly the thermal equilibrium one: 
Ck = Cl'' = and Vk = = {e^"/''^^ + (the 
Fermi function). 

Finally, once the solution of the kinetic equation ((26l) is 
obtained, the density of excitations and energy produced 
during the quench can be expressed as 

^ = ^E^^ (29) 

A:>0 

^ = i^E^^^'^- (30) 

fc>0 

We conclude this section by commenting on a useful ap- 
proximation to evaluate numerically the kernel of Dqk, 
Eq. (|27|) . discussed in full detail in Appendix A2. It 
consists in approximating the evolution operator Uk ap- 
pearing in Dqk with Ukit,t — s) ~ cxp (^i'Hkit)s^ , thus 
obtaining 

Dqk ~ i J dsg>{s) cxp (~i'Hk{t)s^ cxp (i'Hk{t)s^ . 

(31) 

This is again consistent with the separation of time scales 
mentioned above, and in particular with the Markov ap- 



where g[E] is the Laplace transform of bath correlator 
Gk^ is the thermal equilibrium value of the pop- 



proximation. Indeed, while the exact relaxation rate ma- 
trix (P7)) depends on the velocity of the quench, if the 
quench is slow on the time scale characteristic of the bath, 
the correlation fimction g''{s) can be seen as strongly 
peaked at s = 0. Hence the system can be considered 
"frozen" at the instantaneous value of h{t) and, consis- 
tently, its evolution operator is the exponential of the 
Hamiltonian. 



V. RELAXATION TIME 

In order to make further progress in understanding the 
quench dynamics of the system we will first extract from 
the kinetic equation the characteristic relaxation time for 
the populations of the excitations Vk (see Eg. ipS]) ') as 
a function of the magnetic field h and the temperature 
T. For this purpose, it is sufhcient to consider only the 
diagonal elements of Eq. ([28|) . This is equivalent to the so 
called "secular approximation" for the master equation^, 
which is valid for weak couplings (a ^ 1 in the present 
case). For generic N, we deal with a set of equations 
(only A^/2 modes are independent) of the form: 

jVk =ak+Y^ bkqVq + ^ CkqPkVq . 

1 q 

The asymptotic relaxation can be studied by linearizing 
the previous set of equations near the thermal equilib- 
rium fixed point. We obtain, with the vector notation 
SV = {Sri,SV2, . . . , SVN/2T' where SVk = Vk - T'f : 

^Sr = -n 5V (32) 

where non-linear terms in SVu have been neglected. The 
diagonal and off-diagonal elements of the N/2 x A^/2 ma- 
trix TZ are: 



(33) 



(34) 



ulation of the ground-state of mode fc, Q*j!^ = 1 — P^'*, 



I 

7^fefc = ^ E [gf{^-cos{Bk + eq))g[-Ak-^q\+gf{l + cos{ek+eq))g[Kk-^q] 

q>0,q^k 

(1 - cos(0fc + Oq)) g[Kk + A,] + (1 + cos(0fc + dq)) g[-Kk + A,]] 
+ I 4sin2 9k (t/f <?[~2Afe] + T'f 5[2Afe]) 

Hkq = ^ (1 + co^iOk + Oq)) g[~Kk + A,] + gf (1 - cos(0fc + Oq)) g[~Kk ~ A,] 

-T'f (1 + cos(0fc + Oq)) g[Kk - A,] + (1 - cos(0fc + Oq)) g[Kk + A,]] 

I 
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and 



arccos - 



(cos k + h) 



(35) 



The eigenvalues of TZ, {Xi}, are the characteristic relax- 
ation rates of the long-time dynamics. Hence the solution 
of Eq. ([5^ for each population would be a linear combi- 
nation containing all the characteristic relaxation times: 



Tkjt 



At long times t ^ (minj Xj) t all modes relax with 
the same relaxation time r. In the following we first an- 
alyze the longest relaxation time r, extending the results 
presented in Ref. [s^; we then study the structure of the 
entire spectrum of relaxation times. 




Figure 4: Relaxation rate 1/r as a function of T and h for 
iV = 400 (here 7 = l,and s = 1). 



In Fig. m we show the general behavior of r in the 
finite temperature phase-diagram, calculated by numeri- 
cally diagonalizing the matrix TZ. As T ^ 0, r diverges 
and close to the critical point two different behaviors are 
found in the semiclassical regions and in the quantum 
critical region (see also Fig[T]): 



T ^ (X 



(36) 



These relations extend the results obtained for the relax- 
ation rate in Ref. [s^ to the generic case of non-Ohmic 
baths and give the exponent 9 as 9 = 1 + s. 

An analytic expression for the power-law scaling inside 
the quantum critical region can be obtained by approx- 
imating the smallest eigenvalue of TZ with the smallest 
diagonal element. This is justified by the fact that the off- 
diagonal elements are of the order 0{1/N) (see Eq. (|34p ). 
For h = he = 1, considering the gapless mode k = n, we 



have from Eq. (|33p in the continuum limit: 

- (c^fffhA^l+PfgiA,]) 
Jo 

^"X '^inh(A,/r) 



diag 



8a(l-2-i-^)r(l + s)C(l + .s) (7/T)- 



(37) 



where F and Q are the Gamma and the zeta functions, 
and we used the critical dispersion relation A^ ~ 7(7r — 
q) obtained by linearizing Eq. ([T6|) around the gapless 
point k = TT (we extended the integration to —00 since at 
low-temperature only the low-energy modes contribute 
to the integral). Fig. O demonstrates that the analytical 
expression in Eq. (|37p agrees very well with the numerical 
solution (obtained by diagonalizing TZ), especially at low 
temperature. 
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Figure 5: Relaxation rate 1/t as a function of T as obtained 
from the exact diagonalization of TL (symbols) and the ap- 
proximation in Eq. (|37p (dashed lines). Upper panel; s = 1, 
3, 6, with 7 = 1; lower panel: s = 1 with different 7 = 0.3, 
0.5, 1 from top to bottom. 

As we have shown in Fig. [5l inside the quantum crit- 
ical region the exponent 9 is universal within the range 
of anisotropy < 7 < 1 where the system belongs to the 
Ising universality class. This suggests a relation between 
9 and the critical indexes of the quantum phase transi- 
tion. Indeed, it can be shown, within the Fermi golden 
rule (sec Appendix B), that for a generic system coupled 
to a bosonic bath the following expression holds inside 
the quantum critical region: 



1 ^ rj^S^djz 



(38) 
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An important feature of the relaxation dynamics can 
be extracted by analyzing the spectrum of the eigenval- 
ues {Xj} of TZ. In Fig. [6] the Aj's are shown for some 
values of temperature and magnetic field. We find that 
in the semiclassical regions T ^ A the smallest eigen- 
value of TZ is separated from the rest of the spectrum by 
a gap (even in the N oo limit). On the contrary, in- 
side the quantum critical region such eigenvalue merges 
with the rest of the spectrum. This can be quantified by 
the relative gap of the spectrum of relaxation times, that 
is identified by (Aj^^ — Aj~^)/Aj~^, being A1.2 the lowest 
eigenvalues of TZ (see Fig. |V|. Such result indicates that 
while the exponential divergence of the relaxation time 
T cx exp{A/r} in the semiclassical regions is due to an 
isolated eigenvalue, the long-time behavior in the quan- 
tum critical region is, instead, built up by a continuum 
of eigenvalues contributing to the r^^ cx scaling. 
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Figure 6: Spectrum of the eigenvalues of TZ for the Ising 
model (7 = 1) with Ohmic baths (s = 1), for iV = 100. The 
values of temperature T = 0.1 and magnetic field h = 0.8, or 
1 are chosen to belong, respectively, to the semiclassical and 
the quantum critical region. 



VI. ADIABATIC QUENCHES 

Equipped with the kinetic equation and the knowledge 
of the scaling of the relaxation times, we now analyze the 
quench dynamics of the model in Eq. (jlSp ) by solving the 
kinetic equation (|26p numerically. The system is initial- 
ized at hi ^ he in equilibrium with the bath at a fixed 
temperature T <^ A(ft.i), and the transverse field is then 
ramped linearly h{t) = hi — vt down to a final value hf 
(the bath temperature is kept fixed). 

In Fig. [S]we plot the density of excitations as a func- 
tion of the quench velocity for different system sizes (a 
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Figure 7: Ising model (7 = 1) with Ohmic baths (s = 
1). Relative gap between the two longest relaxation times: 
(A2 ^ — Aj^^)/Aj^^ (Ai,2 being the two lowest eigenvalue of TZ); 
crossover lines T — Tcross = 1| are plotted for comparison. 
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Figure 8: Density of excitation £ (circles) Vs quench velocity 
for different system sizes iV = 26, 50, 100, 200, 400, 800 
from bottom to top, according to the arrow; the points cor- 
responding to A'^ = 800 and A'^ = 400 are indistinguishable. 
Parameters are set a = 0.01, T = 0.1 7 = 1 and s = 1 and 
the quench is halted at ft/ = 0.8. Dotted line is the coherent 
contribution Ecoh obtained for a = and stars represent the 
incoherent contribution £inc due to thermal excitation (see 
text); both curves refere to N = 800, even if for Sine the 
same curve is obtained already at A'" ~ 30. 



similar behavior is obtained for the density of energy). 
Additionally we considered separately the coherent {Ecoh) 
and incoherent contribution {Sine) to the final density of 
excitations. The first one is obtained by integrating the 
kinetic equation for a = 0, i.e. no coupling with the 
bath. The incoherent term is due to thermal excitations 
created by the bath and it is obtained by integrating 
the kinetic equation and ignoring the unitary evolution 
term i\Hk,G^], responsible for the coherent excitation 
process. 

In order to understand the two excitation mechanisms 
we analyse directly the dynamics of the populations 7^^. 
From the results shown in Fig [2] (left) it emerges that 
excitations are generated close to the critical point and 
when the system is driven in the semiclassical region (and 
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Figure 9: Populations of the excited states Vk for A'^ = 400, 
V — 0.0017 and the same values of Fig. [8] Left: dynamics 
of two low energy modes ki ~ n (dashed line) and k2 ~ 0.97r 
(dotted-dashed line) obtained by solving the kinetic equation 
(instantaneous thermal equilibrium values are plotted as ref- 
erence as dotted lines). The Inset shows the energy levels 
near the critical point and the scale of temperature; marked 
levels of the excited band refer to ki and k2. The energy gap 
closes at = TT. Right: distribution of Vk as a function of 
the mode k at h{t) = 1. Stars represent the excitations cre- 
ated incoherently and triangles are the coherent excitations 
produced in the case of no coupling to the bath. The two 
excitation mechanisms act on different energy scales, being 
lowest energy modes coherently populated and the highest 
one thermally excited. 




Figure 10: Density of energy (lowest panel) and of excitations 
versus quench velocity v for hf = 0.8, and 1. Parameters are 
set to 7 = 0.7, s = 1.5, a = 0.01, and T = 0.15, or 0.1 
(upper and lower curves of each panel). Circles are obtained 
by solving the kinetic equation; dotted lines are the coher- 
ent contributions Scoh and Ecoh evaluated by solving the ki- 
netic equation for q = 0; a fit gives correctly Scoh oc y^ij and 
Ecoh OC V consistently with the KZ scaling-law for excitations 
^ and the modified scaling we derived for the energy den- 
sity (O. Solid and dashed lines are Eqs. ^ and ^ using 
for the incoherent contributions the expressions (|39|) and (|40p 
and their linearized forms (|41|l and (|42p . respectively. 



T ^ A) they are relaxed out by the bath. The density 
of excitation generated is the sum of the the incoherent 
and coherent contribution, thus proving the validity of 
the Ansazt ^ and (g]) (see Fig. [l (right)). 

For the XY model the integration of Eqs. ^ and ^ 
can be performed explicitly by using the critical spectrum 
Afc 7(7r — k). We obtain: 



f. — 



log 2 



2ttj 
247 



(39) 
(40) 



where the latter holds for quenches halted at hf = he- 
In the previous formulas, the expression derived for r in 
Eq. ([57)1 can be used to get a fully analytical expression. 
Expanding the exponentials in Eqs. ((39|) and (|40|) we ob- 
tain 



8 log 2 



(p{s) 



-2-s „.-i 



(41) 
(42) 



where (^(s) = (1 - 2-^-'')r{l + s)C{l + s). The previous 
relations are consistent with Eqs. pH)) and pT|) with 9 = 
1 + s (see Eq. ([38l)). 

In Fig. [TO] the density of excitations and energy ob- 
tained from the solution of kinetic equation is compared 
with the scaling-law derived in Sec. [Til using the specific 
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Figure 11: System and bath parameters are fixed as in of Fig. 
1101 Upper panel: Data collapse of fi„c and Ei^c obtained 
from the kinetic equation; data refer to 10"'^ v < 10~^ 
(data relative to hf = 1 for Sine are rescaled by a factor 2, 
since in this case only half quantum critical region is crossed). 
Lower panel: scaling of Vcross is obtained equating Sine = Scoh 
and analogously for E. The fits confirm the scaling predicted 
by ([TU)) . (fTTj) and {T^J, (O, that, for the specific case of 
s = 1.5 considered, are shown in their corresponding plots. 
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Figure 12: Density of excitations versus quench velocity v for 
hf = 0.8, 0.6, 0.4, 0.2, 0. Parameters are set to 7 = 1, s = 1, 
a = 0.01 and T — 0.1. Upper solid line is the scaling-law Q 
using the expression (|39p . Decreasing the value of hf, the 
crossing of the semiclassical region, after the critical point, 
becomes more relevant at low v and scaling no longer holds 
strictly. 



expressions, Eqs. ([55)) and derived above for the 

XY model. The scaling-laws are found in good agree- 
ment with the numerical data. The results shown in Fig. 
111! further confirm the scaling as a function of the tem- 
perature and the relations p^ . for the crossover 
velocity. 

Finally we comment on the role of the final value of 
magnetic field hf at which the quench is halted. The 
agreement with the scaling Ansatz becomes worse for de- 
creasing hf (sec Fig. [T^. This is due to the non-critical 
relaxation induced by the bath when the system crosses 
the semiclassical region after the critical point (see Fig.[T] 
and Fig. [S] left). At low v the time spent therein at a rel- 
ative low-temperature T ^ A is so long that the bath is 
able to relax the excitations created close to the critical 
point. 



To check the validity of the scaling-laws, we considered 
the specific case of the quantum XY model pi]) coupled 
locally to a set of bosonic baths, Eq. (|17p (see Fig. [2]). In 
order to study the dynamics we derived a kinetic equa- 
tion, within the Keldysh formalism. A detailed analysis 
of the characteristic relaxation time obtained from the 
kinetic equation was given in Sec. [V] An analytic ex- 
pression for the critical relaxation time was obtained in 
Eq. (j37[) and verified in Fig. [5l As shown in Appendix 
B, the scaling of the latter as a function of the tempera- 
ture is related to the critical exponents of the model (see 
Eq. (j56p ). Finally, we considered the quench dynamics. 
The kinetic equations derived allow us to study the dis- 
sipative dynamics also beyond the universal regime. We 
checked the scaling-laws derived and their range of valid- 
ity in Figs. [TO] and [TH 

We remark that the method described here to obtain a 
kinetic equation for the XY model, may be extended to 
describe the dissipative dynamics of other models that 
can be mapped into fermionic degrees of freedom, like 
other spin chains and ladders or certain 2d models of the 
Kitaev-type. 
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APPENDIX Al: KINETIC EQUATION 



VII. CONCLUSIONS 

We have studied the dynamics of a quantum critical 
system coupled to a thermal reservoir and subject to an 
adiabatic quench across its quantum critical point. We 
considered the regime of weak coupling, low-temperature 
and slow quench velocity. 

The bath has two effects on the system: the first one 
is to create excitations inside the quantum critical region 
and the second one is to trigger the relaxation of the exci- 
tations created close to the critical point when the system 
is driven in the semiclassical region (see Fig. [T]). While 
the first mechanism is universal, being entirely ruled by 
the critical properties of the low-energy spectrum, the 
latter depends on the details of the system far-off the 
critical point. Hence, as far as the evolution is halted 
close to the critical point and the non-critical relaxation 
mechanism is negligible, universal scaling behavior is re- 
covered. We derived scaling-laws for the density of energy 
produced by the quench at finite temperature extending 
the previous results obtained for the density of excita- 
tions in Ref. [s^. 



Here we present a detailed derivation of the kinetic 
equation. Apart from the Keldysh pT|) . lesser (|23p and 
greater ([24]) Green's function we need also the retarded 
and the advanced ones, and also the bath Green func- 
tions: 

9^^^Htl,t2) = (-)z(X,(i2(l))X,(ti(2))) 

9f'^\h,t2) = (-)*e(t2(i)-ii(2))([^,(<i), X.ih)]) 

where we used commutators (anticommutators) for the 
retarded and advanced bath (system) Green's functions. 
The starting point is the Dyson's equation p2[) . In 
order to obtain from the Dyson's equation ([^^ an 
equation for the lesser and greater Green's function we 
use the Keldysh book-keeping for a generic convolution 
C{h,t2) = J^dtA{h,t)B{t,t2) is C'-('')(ti,i2) = 

/pdiA'-(")(ti,i)B'-(")(i,t2) and C<>\h,t2) 

J^dhA-{h,t)B<(>^{i,t2) + A<i»{h,i)B-{t,t2)^. 
Using the previous formulas we rewrite the Dyson's 



11 



equations as: 

^at,G<(^l,^2) = nkih)Gk{ti,t2) + (43) 
/ dt ^l(ti , t)G< {t, t2) + J:< ih , t)Gl it, t2) 

Jo 

and an analogous equation for dt2- Wc are interested in 
the equal-time Green's function and hence we perform a 
change of variables: 



t = 



tl+t2 



St = ti- t2 

whose Jacobian is simply dt = 9ti + and dst = 
^{dti —dt^)- At equal time {5t ~ 0), for the lesser Green's 
function = G<{ti,ti) = G^(t) we get: 

idtG< = [Ukit], Gk] + (44) 

^k''^k + ^fe ■ ^ ^k ' ^k ~ '-'k ' ^k 

where the dot indicates the convolution SJJ • = 
J* dt S5J(t, t)G^ (t, t). Now wc use the relations: 

G'\t,,t2) = 0{t,-t2){G>{tut2)-G<{t,,t2)) 
G"(tl,t2) = 0(t2-ti){G<iti,t2)^G>{ti,t2)) 

and similar relations that hold also for the S*"'" (see^): 

S^('')(tl,<2) = S'(5(ti,t2) + 

e(ii(2) - ^2(1)) (s>(<Hii,t2) - S<(»(ti,i2) 

where we can neglect the term E*^ that only rcnormal- 
izes the Hamiltonian and is not relevant in our case (see 
Eq. ([^5]) below). At equal times we get: 

za,G< = [7i,,G<] + 

v> v< -L n< v> v< 

k ' ^k fe ' ^k ^ ""fe ' fe ^k ' k 

We now perform a Markov approximation. This will 
transform the integro-differential kinetic equation into a 
differential equation. Let us define the interaction picture 
for a general function: 

dk{tiM) = ul{ti)Ok{ti,t2)Uk{t2) 

where lAk is the free evolution matrix for the system obey- 
ing iUk = TikUk- Such transformation gauges away the 
free evolution and the new Green's function Gk dynamics 
is solely governed by the self energy: 

zdtG< ^ S> • G< - E< • G> + G< • E> - G> • E< . 



Since the self-energy carries the "small" perturbative 
coupling parameter the evolution of Gk can be regarded 
as "slow" with respect to the time scale of the self energy, 
that is the same as that of the bath. In fact S contains 
the bath-correlation function g{ti,t2) (see Eq. p8)) ) that 
is strongly peaked at ti ~ t2 because of the assumption of 
a cutoff-time for the bosonic modes (see Sec. lIVp . Thus 
we can take G out of the convolutions: 



lotu,: ^ ( f'^dit>{t,t)^ G< - (^jytt<(t,t) ) Gl 
dlY> {i, t)^ - G> (^J^ dtt< {t, t) ] (45) 



Eq. (|45|) is quite general and it is based solely on the as- 
sumption of Markovian baths. Wc now use the explicit 
form of the self-energy for the coupling system-bath ([^0]) 
with / = 1, that within the self-consistent Born approxi- 
mation reads: 



Sfc(ti,t2) = gfc-,(ti,t2)r^G,(ii,t2)r^ 

9 

= -^g(ti,t2)r^E^9(^i'^2)r^ (46) 



where gq{ti,t2) — ~i {T Xq{ti)Xq{t2)) is the non- 
interacting bath Keldysh Green's function that does not 
explicitly depend on the moment q (since all baths have 
the same spectral function) . In Eq. we neglected the 
polaronic shift contribution (corresponding to the tad- 
pole diagram, Fig. [3}j) 



^i{ti,t2) ^ -^S{tit2)T' f dtg{t^,t)Y,Ty[T''Gq{lt)] 

-'-y q 

(47) 

In fact, being such term proportional to a (5(ii.t2), it 
has only the irrelevant effect of renormalizing the Hamil- 
tonian (see Sec. IIV|) . Using again the Keldysh book- 
keeping^Sii^, we obtain from Eq. pS)) . for the lesser and 
greater self energy 



Sf(ii,t2) = ^3^(ti,<2)r^E^f(*i'^2)r^ (48) 

9 



(notice that g^{ti,t2) = — .g^(ti, i2)*)- Evaluating ex- 
plicitly the self-energy kernels we obtain: 
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N 



/ dt g> (t - t)Ul {t)r'U, {i)G> {t, t)Ul it)r'l^k {t) 



^ J2i^l{t)T^U,{t)G>{t,t) Q^°°ds g> {split - s)T^Uk{t - s) 



N 



J2 / dig>{i,t)ul{T)T^G>{i,t)T^Uk{t) 



N 



J2 / dt g> (i - t)ul {T)t'U, it)G> (t, t)Ul (tyUk (t) 



^J2i L 9^i-')Mit - ^y^l^iit - ) G>it,t)Ulit)r^llkit) 



r 



where ~ refers again to the Markov approximation. For 
the greater kernels simply interchange "<" with '>". 
Finally, in the Schrodinger picture, using the relation 
G> = -il + G<, we obtain Eq. 



APPENDIX A2: APPROXIMATION FOR THE 
KINETIC EQUATION MATRICES D 

In this appendix we comment on the validity of the 
approximation (|3ip for the matrices (|27p appearing in 
the kinetic equation. To calculate D exactly, we need 
to know the evolution operator Uk, solution of the dif- 
ferential equation iUk ~ UkUk- This can be obtained 
exactly by mapping the dynamics of a generic mode k 
into a Landau-Zener two-level system dynamicsi^*^ 



n 



LZ 



LZ. 



(49) 



with A^'^ = 7 sin A; and hj:^ = —vt that can be obtained 
from (|15|) by a simple change of the variable t. Using 
the solution of the Landau-Zener problem for a quench 
that starts at i = — oo we have for the matrix elements 
oiU^^{-oo,t), 



U^^{~oo,t) 
the following results: 

^n{t) = e*^ exp J - 



U^i^{t) Ut,^{tr 



TT A 



LZ\ 



U^fit) 



A 



LZ 



8v 



V t] 



where p = —i (A^^) /2v and Vp are parabolic cylinder 
functions. Finally the evolution operator from generic t 
to t can be obtained using the simple property: 

U^^{t,t) = U^^{-~oo,t)U^^{i,~oo) 
= U^^ {-oo,t)U^^\-oo,t) 

The second ingredient we need in order to calculate D 
is the bath thermal equilibrium correlation function g^ . 
From its definition: 



5>W 



-i{X{t)X\Q)) 



2 I ^ — iujfjt 



duj J(w)[e-''^*(1 4- neiuj)) + e*"*ns(w) 



where ub = l/(e'^/-^ — 1) is the Bose function and we 
used the definition (fT9|) of spectral function for the bath 
J{lo) = X)/3-^^'^('^ ~ ^p)- The correlation function can 
be written explicitly as: 

/•oo , ^ 

g^{t) = ~i 1 dti' J(ti-') I cotli( — ) cos(a;r) — i sin(ti;T) 

= -2iaT'+^T{s + l) X 



Cis + 1, T ^ ) + C(s + 1, T ^ 



where F is the Gamma function and C(z, u) = 

E:r=o(;iT^-«^o,-i,-2.... 

We are now able to calculate explicitly the matrix D 
and check the validity of the approximation in pip . As 



/2v 



exp ■ 



p 1 ((^1 + *)v^ ^^^^'^'-^ ^cc. IIVI the approximation consists in consid- 
I ering the instantaneous transition rates induced by the 
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bath, independent on the velocity of the quench. This 
is ultimately justified by the assumption of "fast" and 
memory less Markovian baths. Hence, for slow quenches 
when the typical time of the quench is much larger than 
the typical bath time-scale we expect that the magnetic 
field can be regarded as not evolving for the bath. Within 
the approximation in Eq. (j31|) we can perform explicitly 
the integration for the matrix elements of Dqk over time, 
giving the Laplace transform of the bath Green function. 
We are interested only in the real part of the latter 



since the imaginary part gives a renormalization con- 
tribution that is negligible in the weak coupling limit 
a — !■ 0^. In the basis of the eigenvectors of Hk we ob- 
tain: 



g[E] - *.g>(t)e'^^*] 
= 7t{J{E)-J{-E)) 



exp iPE) 
exp iPE) - 1 



(50) 



jjapprox^}_f cos++.g+ + COS+ g++ i (sin+ g + + sin++ .g )\ , - 

9fc ^ 2 V -i (sin+" + sin++ cos~+5 -cos++,g"+ J ' ^ ' 



where we defined 

cos^* = ± cos ± cos 9q 
sin^* = ± sin 6'fe ± sin 0g 
5±± = g[±Afe±A,] 

The sin and cos are geometric factors specific of the sys- 
tem operator that couples to the bath (in our case cr^), 
while the Laplace transform of the bath g[.] (see [50)) car- 
ries information about the relaxation rates between the 
different energy levels, and depends explicitly on the tem- 
perature and on the nature of the baths (i.e., its spectral 
function). For simplicity we consider the equal indexes 
Dkk matrices. The same results hold also for unequal in- 
dexes since the integral of the matrix elements have the 
same structure in both cases. In Fig. IVIII we compare 
the matrix elements obtained using the exact evolution 
operator with the ones given by Eq. ([5T|) . The agree- 
ment is good, validating the approximation. Deviations 
appear only in the limit of fast quenches -y/w 3> T, and 
in such regime the bath has a less relevant effect on the 
dynamics because of the short interaction time during 
the quench. Besides that, deviations appear far from the 
critical point (corresponding to h'"^ ~ 0), i.e., far from 
the most relevant part of the quench according to Secs.HIl 
andlVD 



APPENDIX B: FERMI GOLDEN RULE FOR THE 
RELAXATION TIME 

In this section we derive an expression for the crit- 
ical relaxation time using the Fermi golden rule for a 
generic system interacting with a bosonic bath. Let us 



assume the system-bath interaction Hamiltonian to have 
the form Hint = where A and Z are system and bath 
operator respectively. Consider a quench of the system 
from zero temperature to a certain finite T. The transi- 
tion rate for the process of thermalization in presence of 
the reservoir pf {\GS) {GS\)s ^ ® pf (where B 
and S refer to system and bath, respectively) is: 

- = 2^Y.^^^f+^k-E.,-EGs)PB{E^/T) 

xP's\Ek/T)\{k, f\H,nt\GS, (52) 

where i, / and k address the bath eigenvalues and 
the final state of the system respectively; -Ps{b) 
are thermal weights. We rewrite the (5-function as 
^ r°° dte-^(^/--^')*e-'(^''-^Gs)t. Summing over / and 

i we get the bath correlation function z{t) ~ {Z{t)Z{0)): 
- = y / dte-*(^^-^'5^)*z(i)P|''(-Efc/r)| (fc| A \GS) |2 

(53) 

The time integral gives the Fourier transform of the bath 
correlation function, that we parametrize as 

z[E] = J{E)f{E/T) . (54) 

For instance, for a bosonic bath with spectral function 
J{E) cx E'^ we have 

{j{E){l + nB{E/T)) E>0 
^ \ji\E\)nBi\E\/T) E<0- 

By integrating over the k-modes (setting Eqs = 0)i using 
the critical density of states p{E) oc E'^^^~^, we get: 

i (X jdEp{E)z[E]Pf{ElT)\AGs{E)W 
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Figure 13: Matrix elements of Dkk as a function of the 
rescaled field h^'^ (|49p : ft^'^ = correspond to the critical 
point /i ~ 1 for the relevant low energy modes. Lower panel; 
diagonal elements iiDkk]^! (^P) ^'^'^ ~*(^fcfe)i2 (down); up- 
per panel shows the difference (Phk)x\ ~ (^*:fe)22- Plots refer 
to T/A^^ = 5; continuous line is the approximation (|5ip 
(which is independent on v) and symbols are the exact value 
for y'Ti/T — 0.1, 1, 5; deviations from the approximation are 
appreciable only for the last value of v. 



where Acs{E) = {k{E)\A\GS). Now, assuming that 
the low-energy modes (that are the relevant ones at 
low-tcmperature) arc coupled uniformly by the bath 
AcsiE) ~ AasiO) we obtain 



-oc dEE'^/'--^ J{E) f{E/T)Pf{E/T) (55) 



and finally, using for the spectral density J{E) oc E^ and 
performing a change of variable to x = E/T we obtain 



cx T^+'^Z^ . (56) 
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